Thermal conductivity in superconducting borocarbides LuNi 2 B 2 C and YNi 2 B 2 C0 
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We have recently proposed the s+g wave model for su- 
perconducting borocarbides. In spite of a substantial s-wave 
component, this order parameter exhibits the y~H dependent 
specific heat and a thermal conductivity linear in H in the 
vortex state. This is characteristic for nodal superconductors 
when T,r <C A where F is the quasiparticle scattering rate 
and A the maximum superconducting gap. Here we investi- 
gate the thermal conductivity parallel to the c- and a- axis 
in a magnetic field tilted by 9 from the c- axis and rotating 
within the a-b plane. 
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The superconductivity in the rare earth boroc 
bides LUM2B2C and YN12B2C is of great interest!! 
We have proposed recently the superconducting order 
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where i? and ip are polar and azimuthal angle of k re-, 
spectively. Recent thermal conductivity experiments^! 
suggest that crystallographic [100] and [0101 are the nodal 
directions, i.e. the order parameter of Eq. in is rotated by 
? in the a-b plane. This gap function accounts for the 

\f~H~ dependence of the specific heat and the H- 1U 
term in the thermal conductivity observed recently 
The aim of this paper is to generalize an earlier result^ for 
Kzz, Kxx and K xy for general magnetic field (H) orienta- 
tion given by the polar angle 6 with respect to the c- axis 
and the azimuthal angle cj). First in the absence of H the 
specific heat and the electronic thermal conductivity for 
r < T < A are given by 
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where Co=(^j-) 2 [ln(2 y ^)] 2. Note that n xx obeys 

the universal behaviour while k zz does not. This is 
because the heat current operator j z vanishes on the 
four second order nodal points (■#, ip) — (5,±f) and 
(f ,±3f ) for A(k) given in Eq.(|). Also this leads to 

arlt ln( . f\ r ) dependence of k zz as discussed below. 
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FIG. 1. Angular dependence of I zz (8,(j}) which determines 
Kzz(0,<j>) in the superclean limit for s+g wave and d- wave 
case (up to log- terms in Eq. 3). Note the different scale in 
the two cases. 



In the presence of a magnetic field with general ori- 
entation defined by (0, <j)) the specific heat and-thermal 
conductivities in the vortex phase are given byEnlj 
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Here we have assumed the superclean limit defined 
by V Ar <C vV eH with v = ^/v a v c where v a , c denote 
the anisotropic Fermi velocities. The angular depen- 
dences of k zz and K xy according to Eqs. ElQ for the s+g 
wave case are shown in the left panel of Figjl] and in 
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FIG. 2. upper panel: polar angle variation of I zz (0,45) for 
s+g- and d- wave case, lower panel: Angular dependence of 
the product I-(8,<f))- l+(0,(f>) for s+g case which determines 
the thermal Hall coefficient K xv (9,cj>). It vanishes for field 
halfway between the nodal directions due to current compen- 
sation. 



Fig.^j. For comparison, we also present the correspond- 
ing angular dependence of n zz for the d^— 2 state wi 



A(k) =Acos(2</>) as irj-4wth T c cupratesH, CeCoIn 5 ti 
and K-(ET) 2 Cu(NCS) 2 liM. Of course for d- wave su- 
perconductors the universal zero- field behaviour is valid 
both for k xx and for k zz , and both exhibit a similar an- 
gular dependence in the vortex phasdI3. In this case the 
dependence on field angles 0, </> is given by 
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cusp like feature develops for 0=90° and <f> = ± 45° due 
to the (second order) point node, while in the d x 2_ y 2 
wave case with an extended line node along c no cusps 
appear and also the absolute value of angular variation 
is much smaller. This is clearly visible from the upper 
panel of Fig. |^ which also shows monotonic 0- dependence 
for s+g- wave and nonmonotonic behaviour for d-wave. 
The latter has a minimum at 9 m ~ 47° which is due to 
a maximum Doppler shift for 0=45° resulting in a dom- 
inating term I zz (6,(/>) ~ 1 - (5/64) sin 2 (20)+.. . Note 
that l zz (9, 4>) in Fig. [I] exhibits a rather sharp minimum 
as function of 4> at 9 m wheras for 9= 90° the minimum 
is flat. Experimentally however the n zz thermal conduc- 
tivity shows very strong cusps at 0=902 (and 4>—0 due to 
rotated order parameter) in YNi2B2CD. This is a strong 
point for the s+g wave case being the appropriate one for 
YNI2B2C and LUM2B2C . Therefore the thermal conduc- 
tivity in the superclean limit can discriminate s+g- wave 
against d- wave superconductivity. 

The angular dependence of the thermal Hall coefficient 
K xy in the s+g wave case is shown in the lower panel of 
Fig. H It exhibits a sign change as function of and 
varies smoothly with 9. In the d- wave case n xy looks 
rather similar. ■— ■ 
As we have already discussed elsewhereEJ the thermal 
conductivity provides a unique window to look at the 
nodal structure of A(k) in unconventional superconduc- 
tors. 
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Then k tt and 



Eq. ^ but now for d- wave: 
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The (f>- dependence of n zz for various is shown in 
comparison to the s+g wave case in Fig. |IJ As is read- 
ily seen from Fig. |l| in the s+g- wave case a pronounced 
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